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DEVELOPMENT OF SONIC -BOOM SIGNATURES IN A 
STRATIFIED ATMOSPHERE 

By Raymond L. Barger 
Langley Research Center 

SUMMARY 

Equations describing the development of the pressure signature of a supersonic pro- 
jectile or equivalent body of an airplane in straight steady level flight in an atmosphere 
having vertical gradients of temperature, pressure, and wind have been derived. The 
analytical method utilizes results of ray-tube calculations for a stratified atmosphere in 
the context of an analysis by Whitham. 

INTRODUCTION 

The sonic-boom literature contains many treatments, which use the acoustic weak 
shock approximation, of the effects of nonuniform atmosphere on the boom intensity. 

(See, for example, ref. 1.) However, this acoustic approximation treats the wave as a 
single pulse and does not account for any details of the nature or development of the pres- 
sure signature. Even in the uniform atmosphere limit this approximation predicts an 
incorrect variation of amplitude with radial distance from the flight path. An attempt to 
apply this linear acoustic theory to the problem of calculating the pressure signature of a 
supersonic body of revolution in an atmosphere with a temperature gradient has been 
described in reference 2. However, this linear theory fails to predict the formation of 
shock waves and other nonlinear features of the signature development. 

On the other hand, these nonlinear phenomena are predicted by the theory (refs. 3 
and 4) that Whitham has developed and applied to the calculation of signatures in a uniform 
atmosphere. The present analysis represents an extension or application of the nonlinear 
theory of Whitham used with a ray-tube calculation procedure similar to that used in ref- 
erence 5 to describe the development of the signature in a stratified atmosphere with 
horizontal winds but without vertical winds. 

SYMBOLS 

The symbols a, p, u, v, X, p, v, and p without a subscript or tilde apply 
for conditions undisturbed by the wave, z units below the flight level. 
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defined by equation (14) 
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Mach number 

unit vector normal to element of wave front associated with ray under 
consideration 

local pressure 

distance along ray path 

source distribution due to supersonic body 

time at which ray under consideration intersects horizontal plane z units 
below flight altitude 
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T 

u,v 

U 

V 


time variable used in equation (4) 

wind components relative to X,Y,Z coordinate system 

particle velocity in wave of finite amplitude 

flight velocity vector relative to air at flight altitude 


V = 


V 


w 


wind velocity vector relative to air at flight altitude 


w n = w • n 


X,Y,Z right-hand rectangular coordinate system moving with atmosphere at flight 

altitude, with Z-axis directed vertically downward and X-axis directed 
opposite to direction of flight; origin located at nose of generating body 
at time t = r = 0 


x, y,z coordinates referred to X-, Y-, and Z-axes, respectively 

w* wind velocity vector at flight altitude relative to ground 

u*,v* components of w* referred to H,H system 

0 = \Jm 2 - 1 

y ratio of specific heats 

9 angle of inclination of ray with horizontal 

A,p,y components of n 

!// angle of rotation of H,H coordinate system relative to X,Y system 

p density 

t time at which element of generating body, initially at position x,0,0 crosses 

origin of X,Y,Z system 
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a) angle between vertical plane through flight path and plane determined by flight 

direction and initial ray direction 

E,H coordinate system fixed relative to earth with origin directly below nose of 

generating body at t = r = 0, and with E-axis directed opposite to direction 
of flight relative to ground 

|,?7 coordinates referred to E,H system 

£ dummy integration variable 

Subscript: 

0 undisturbed conditions at flight level 

ANALYSIS 

General Considerations 

In deriving the equations for calculating the pressure signature, the coordinate sys- 
tem used is assumed to have its origin at flight altitude and to be moving with the atmos- 
phere at this level so that in this system there is no wind in the plane Z = 0. The X-axis 
is directed opposite to the direction of flight, and the origin is located at the nose of the 
body at time t = 0. Subsequently, this system will be related to a system fixed relative 
to the ground. 

Throughout the analysis the atmospheric variation is assumed to be slight over a 
distance of the order of the length of the signature. If this condition is not applied, cer- 
tain assumptions of geometrical acoustics fail. However, when the atmosphere is such 
that the condition does not hold, turbulence is likely to occur and give rise to a different 
type of sonic-boom distortion which is not described by the present analysis. 

Basic Equations 

The equation for the rays associated with the acoustic wave fronts can be determined 
by a procedure analogous to that used in reference 3 for propagation in a uniform atmos- 
phere. The basic concept is to start with the equation for a sound ray as determined by 
the assumption of ordinary acoustic propagation, and then modify that equation by replacing 
the expression for the undisturbed sound speed by the actual speed of propagation as 
influenced by the finite overpressure in the wave. For any pressure signature there is at 
least one point of zero overpressure, to which corresponds a ray associated with a wavelet 
moving at the undisturbed sound speed. Then this ray is subject to the ordinary acoustic 
refraction laws, as described, for example, in reference 6. The following approximation 
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procedure is used: The ray associated with a point of zero overpressure is adopted as a 
"typical" ray for the signature for the purpose of determining the wave-front normal 
direction and the ray-tube cross-sectional area. Then the nonlinear effect of the motion 
of the other points of the signature relative to the point associated with this ray is deter- 
mined by adjusting the speed of those points to account for the finite overpressure asso- 
ciated with them. 

A fundamental expression required in this procedure is the time at which the wave- 
let associated with the typical ray arrives at a level z units below the plane z = 0 



(See eq. (47) of ref. 6.) 


( 1 ) 


In order to account for the finite overpressure associated with some other ray, 
equation (1) is again used, but with the undisturbed sound speed a replaced by the actual 
speed of the wavelet 


(See ref. 3.) 
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The variation due to refraction of the overpressure Ap = p - p along a ray tube has 
been determined in reference 7 (eq. (28) in section HI of ref. 7). With the use of substitu- 
tions from equations (5) and (13) in section in of reference 7, Ap can be expressed in 
the form 
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The corresponding uniform atmosphere expression is 
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(ref. 3, eqs. (7) and (54), where in eq. (54) y should be in the numerator). 
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By comparing the expression for Ap with the previous Blokhintzev expression in 
the limit of a uniform atmosphere, the parameter c is determined to be 

ill/2 

F(r) 


c = y 


P0 a 0 M5 


2(M 2 - 1) 


Therefore, the expression for Ap becomes 


Ap = y 
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( 2 ) 


Here, the factor F(r) accounts for the variation due to the finite amplitude of the 
wave, while the final factor accounts for the effects of refraction. Replacing a in equa- 
tion (1) with a + U and using equation (2) yields approximately 
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The new time variable T defined by 

T = t • 


F(r) f 




0 ^Jpa(a + w n )|w + an JaJ 1 ^ 


(3) 


r z dt; 

J 0 va 


is displaced from the original time variable by an amount that is constant for fixed z; 
that is, T is independent of r. In terms of T, equation (3) becomes 


T = r- cF(t)J Z — ^ /2 (4) 

0 i/jpa(a + w n )|w + an |aJ / 

Derivation of Ray-Tube Area 

Equations (2) and (4) are the fundamental equations for calculating the signature, but 
before they can be used, A and v must be expressed as functions of z. The following 
procedure for determining the ray-tube cross-sectional area A (which is proportional to 
A) is based on the method used in reference 5. 

The rays considered are those associated with the propagation of the Mach wavelets. 
These wavelets propagate as acoustic waves and near the X-axis have the shape of a cone, 
which forms an angle with the axis equal to the complement of the Mach angle. In order to 
specify a particular ray, it is convenient to consider the initial incremental section of the 
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ray cone (fig. 1). Then any given ray can be specified by the angle co that the projection 
of the ray onto the cross section of the cone makes with the vertical. 


The ray tube to be considered is defined by four bounding rays : the ray emitted at 
x = 0 and specified by the angle co, the ray emitted at x = 0 and specified by the angle 
co + Aco, and the two corresponding rays emitted at Ax = V At. (See fig. 2.) For a ray 
associated with the angle co, the constant "trace" velocities Cj and C 2 in the strati- 
fied atmosphere are defined by (ref. 6, eqs. (41) and (42), where uq = vq = 0) 


since 



or, since 


Cl =^a 0 M 


_ ^o_ _ a o 
2 ' “G ' (l - V - *o 2 


/3 cos co 

"0 = m“" 


(5) 


(ref. 5, eq. (14)) then 


Consequently 


aQM esc c o 

C2= l 


c i 

7T- = & Sin CO 

C2 


( 6 ) 

(7) 


Since the vertical component of the wind velocity is assumed to be negligible, the equations 
of refraction are (ref. 6, eqs. (41) and (42)) 


a + Xu + pv = ACi 


( 8 ) 


a + Au + pv = pC 2 


From equations (8) and (9) 



Equations (8) to (10) can be solved for A and p, with the results 


(9) 
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and 
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Then v is determined by 




v 2 - /x 2 


or with a substitution from equation (10) 


v = 1 - 1 + 


c 2 s 


Define the function f(z,w) by 


x „ / n v(z) + a(z)M(C 1 ,C 2 ,z) 

f z,Ci(w),C 2 (w) = ; — 

a(z)y(Ci,C 2 ,z) 


(13) 


(14) 


Then the y coordinate of the intersection of a ray specified by the angle c o with a hori- 
zontal plane z units below the source of emission is (ref. 6, eq. (46)) 




(15) 


The corresponding intersection point of a ray specified by the angle w + Aw is displaced 
in the y-direction by the increment 


dy 

Ay = — Aw 
dw 

dv 

where -p- is determined by an application of Leibniz's rule 
dw 


dw Jp 


(16) 


(17) 


8f 


In order to evaluate this quantity, -L must be expressed as a function of z. The 

ow 

first step is to expand the derivative by the chain rule 

JL = JL 2 a. , §L in 

3 w 3fi 9w dv 8 w 
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The quantities on the right-hand side of equation (18) can be obtained from previous 
relations. From equation (14) 


and 


From equations (5) and (6) 


and 
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In view of equation (21), only 
af 


dfJL 

9C^ 
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dv 

9C 2 


are required to complete the evaluation of 


- — in equation (18). These derivatives are obtained from equations (11) to (13). The 
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resulting expressions are 
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These results may be summarized as follows: When a, u, and v are specified 
as functions of z, then X and p are determined as functions of z by equations (11) 
and (12). In turn, these expressions substituted into equations (19) to (24) determine the 
quantities needed to express -Mj in equation (18) as a function of z. Then the integral 
of equation (17) can be evaluated so that the increment Ay in equation (16) is obtained 
as a function of z. Now consider the ray specified by the angle co, as was the original 
ray, but emitted at a time At units later, so that it is displaced by the distance 
Ax = V At in the x-direction. Inasmuch as the atmosphere is assumed to be stratified, 
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and since both rays have the same initial inclination angle, they are subject to exactly the 
same refractive effects. Consequently, these rays are parallel, and when they intersect 
the plane z units below the flight level, they are still displaced the same distance Ax 
in the x-direction. Consequently, the area of the intersection of the ray tube with this 
horizontal plane is 

Ax Ay = V -pL Ao> At (25) 

oCx) 

The inclination 9 of a ray z units below the flight plane is determined by 

sin 6 = — - — 
ds/dz 


or 


sin 9 = 


av 


|w + an 

Then the actual ray-tube cross-sectional area is determined by 


A = Ax Ay sin 9 


(26) 


with Ax Ay given by equation (25) and sin 9 by equation (26). 

Since for steady motion in a uniform atmosphere, the function A(z) as used by 
Whitham (ref. 3, eq. (52), where X = °°) is defined by 


A = s = z — sec o> 

0 


the coefficient K relating A and A 


is given by 


A = KA 


K = 




/3^Ax Aw 


Consequently 

A = M?. ilZ. sin 6 (27) 

(32 dw 

Equations (17), (26), and (27) define the function A(z) to be used in equations (2) 
and (4) for the signature calculation. This calculation for the pressure as a function of 
time can be performed in a manner similar to that described in reference 8 for the pres- 
sure as a function of distance. The characteristic parameter designated in reference 8 
as a distance measured along the longitudinal axis from the body nose is replaced with r, 
and the quantity x - /3r of reference 8 is replaced with T. 
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Transformation to Ground- Fixed Coordinates 


The y coordinate of the intersection of the ray specified by the angle with a 
horizontal plane can be determined by equations (14) and (15). The x coordinate of the 
intersection point is given by 

x=f u + aX - dC ( 28 ) 

and the time of intersection by 

t = r z a ( 29) 

do av 

Equations (28) and (29) are derived as equations (45) and (47) of reference 6. In 
order to determine the location of the ray-ground intersection with respect to a coordinate 
system fixed relative to the earth (the S,H system), the relationship of this system to 
the original system at the time of intersection must be determined. The E- and H-axes 
are chosen so that the flight velocity vector relative to the ground is directed along the 
negative E -direction. 

Then the S,H origin is located directly below the X,Y origin at t = 0. It is 
assumed that the velocity vector V, which is directed along the negative X-axis and has 
the magnitude of the airplane airspeed, is determined from the airplane’s instruments. 

The wind vector at flight altitude relative to the ground w* is determined from meteo- 
rological sounding data. Then the relation between the X,Y and S,H systems at t = 0 
is as shown in figure 3. At t = 0, the E,H system is rotated relative to the X,Y sys- 
tem by an angle xp such that 

| = x cos xp - y sin xp 

and 

■q = x sin xp + y cos i p 

At the time of intersection, a translation of axes has also occurred because of the motion 
of the atmosphere at flight altitude. Consequently, at time t 


and 


£ - u*t = x cos ip - y sin xp 
rj - v*t = x sin xp - y cos xp 


where u * and 


v* are the components of w* referred to the 


E,H system. 
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CONCLUDING REMARK 


Equations for calculating the pressure signature of an airplane in straight uniform 
flight in a stratified atmosphere have been derived. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., July 31, 1968, 

126-61-06-09-23. 


REFERENCES 

1. Friedman, Manfred P.; Kane, Edward J.; and Sigalla, Armand: Effects of Atmosphere 

and Aircraft Motion on the Location and Intensity of a Sonic Boom. AIAA J., vol. 1, 
no. 6, June 1963, pp. 1327-1335. 

2. Myers, Michael K.; and Friedman, Morton B.: Focussing of Supersonic Disturbances 

Generated by a Slender Body in a Non- Homogeneous Medium. Tech. Rep. No. 40 
(Contract Nonr-266(86)), Dep. Civil Eng. and Eng. Mech., Columbia Univ., June 1966. 

3. Whitham, G. B.: On the Propagation of Weak Shock Waves. J. Fluid Mech., vol. 1, 

pt. 3, Sept. 1956, pp. 290-318. 

4. Whitham, G. B.: The Flow Pattern of a Supersonic Projectile. Commun. Pure Appl. 

Math., vol. V, no. 3, Aug. 1952, pp. 301-348. 

5. Barger, Raymond L.: Some Effects of Flight Path and Atmospheric Variations on the 

Boom Propagated From a Supersonic Aircraft. NASA TR R-191, 1964. 

6. Groves, G. V.: Geometrical Theory of Sound Propagation in the Atmosphere. Jour. 

Atmos. Terrest. Phys., vol. 7, no. 3, Sept. 1955, pp. 113-127. 

7. Blokhintzev, D.: The Propagation of Sound in an Inhomogeneous and Moving Medium I. 

J. Acoust. Soc. Amer., vol. 18, no. 2, Oct. 1946, pp. 322-328. 

8. Middleton, Wilbur D.; and Carlson, Harry W.: A Numerical Method for Calculating 

Near-Field Sonic-Boom Pressure Signatures. NASA TN D-3082, 1965. 


12 





Figure 2.- Ray-tube geometry. 
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